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Abstract 

A search for narrow baryonic resonances decaying into H^vr^ or H^7r+ and their antipar- 
ticles is carried out with the HI detector using deep inelastic scattering events at HERA in 
the range of negative photon four-momentum transfer squared 2 < < 100 GeV^. No 
signal is observed for a new baryonic state in the mass range 1600 — 2300 MeV in either 
the doubly charged or the neutral decay channels. The known baryon H(1530)'' is observed 
through its decay mode into H^7r+. Upper limits are given on the ratio of the production 
rates of new baryonic states, such as the hypothetical pentaquark states H^^^ or Hg^ , relative 
to the S(1530)° baryon state. 
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1 Introduction 



A search for narrow baryonic states in the mass range 1600 — 2300 MeV, decaying into the 
doubly charged S^vr^ and neutral S^7r+ final states and their antiparticle^ll, is presented. The 
search is carried out in deep inelastic ep scattering (DIS) at HERA. 

Various theoretical approaches imSEl based on Quantum Chromodynamics predict the ex- 
istence of exotic baryonic states composed of four valence quarks and an anti-quark, commonly 
known as "pentaquarks". Such states are expected to form a flavour anti-decuplet and are not 
explicitly forbidden within the Standard Model. 

Several experiments have reported evidence for a narrow resonance with a mass around 
1540 MeV decaying into nK~^ and pK^ final states flU. Such a state could be interpreted as an 
exotic strange pentaquark with a minimal quark content of uudds, lying in the apex of the spin 
1/2 (or 3/2) anti-decuplet. On the other hand, a number of other experiments flU, including the 
HI experiment [i5j|, have reported non-observation of the same state. 

Searches for other members of this anti-decuplet are of interest, in particular for the doubly 
strange E'^~ (ddssu) and [uussd) states, which are also manifestly exotic. The NA49 
collaboration [6] reported the observation of two baryonic resonances in fixed targed proton- 
proton collisions at the SPS at the centre of mass energy s/s = 17.2 GeV, with masses of 
1862 ± 2 MeV and 1864 ± 5 MeV, and with widths below the mass resolution of 18 MeV. These 
states can be interpreted as the "E^' {S=-2, Jg = -3/2) and the S°g(S' = -2, J3 = +1/2) 
members of the isospin 3/2 quartet S3/2 in the anti-decuplet. These findings have not been 
confirmed by several other experiments [l7l[8ll9l [T0l[m[T2| l. 

The search presented here is performed using data taken with the HI detector at HERA. The 
particles are identified through their decay into Att^. The established baryon S(1530)° [il3l 
is observed through its decay mode S(1530)° S"7r+. 



2 Experimental Procedure 
2.1 HI Apparatus 

A detailed description of the HI detector can be found in [fT4ll . In the following only those 
detector components important for the present analysis are described. 

HI uses a right handed cartesian coordinate system with the origin at the nominal ep inter- 
action point. The proton beam direction defines the z axis. The polar angle is measured with 
respect to this axis and the pseudorapidity r] is given by 77 = — In tan |. 

The tracks from charged particles are reconstructed in the central tracker, whose main com- 
ponents are two cylindrical drift chambers, the inner and outer central jet chambers (CJCs) [[T5|. 
complemented by a central silicon vertex detector (CST) [16|. The CJCs are mounted concen- 
trically around the beam-line, covering the range of pseudorapidity —1.75 < r] < 1.75 for tracks 

'Unless explicitly mentioned, the charge conjugate states are hereafter always implicitly included. 



4 



coming from the nominal interaction point. The CJCs lie within a homogeneous magnetic field 
of 1.16 T which allows the transverse momentum of charged particles to be measured. Two 
additional drift chambers (CIZ, COZ) complement the CJCs by precisely measuring the z coor- 
dinates of track segments and hence improve the determination of the polar angle. Two cylindri- 
cal multi-wire proportional chambers facilitate triggering on tracks. The transverse momentum 
resolution of the central tracker is a{p^)/p^ ~ 0.005 / GeV © 0.015 [fTTI . Charge misiden- 
tification is negligible for particles originating from the primary vertex and having transverse 
momenta in the range relevant to this analysis. 

The tracking detectors are surrounded by a Liquid Argon calorimeter (LAr) in the forward 
and central region (—1.5 < rj < 3.4) and by a lead-scintillating fibre calorimeter (SpaCal) in 
the backward region JTSl (—4 < rj < —1.4). The SpaCal is optimised for the detection of 
the scattered electroitl in the DIS kinematic range considered here. A planar drift chamber, 
positioned in front of the SpaCal, improves the measurement of the electron polar angle and 
is used to reject neutral particle background from photons. The kinematics of the hadronic 
final state are reconstructed using an algorithm which combines information from the central 
tracker, the SpaCal and the LAr calorimeter. The DIS events studied in this paper are triggered 
by an energy deposition in the SpaCal, complemented by signals in the CJCs and the multi-wire 
proportional chambers in the central tracker. 

The luminosity is determined from the rate of Bethe-Heitler processes ep — > epq, measured 
using a calorimeter located close to the beam pipe at z = —103 m in the backward direction. 

2.2 Selection of DIS Events 

The analysis is carried out using data corresponding to an integrated luminosity of £ = 101 pb^^, 
taken in the years 1996/1997 and 1999/2000. During this time HERA collided electrons 
or positrons at an energy of 27.6 GeV with protons at 820 GeV (1996/1997, 24.8 pb"^) and 
920 GeV (1999/2000, 75.7pb-^). 

The event kinematics are reconstructed from the energy and the polar angle of the scattered 
electron, which is required to be identified in the SpaCal with an energy above 8 GeV. The neg- 
ative four momentum transfer squared of the exchanged virtual photon and the inelasticity y 
are required to lie in the ranges 2 < Q"^ < 100 GeV^ and 0.05 < y < 0.7. The difference be- 
tween the total energy E and the longitudinal component of the total momentum p^, calculated 
from the electron and the hadronic final state, is restricted to 35 < — < 70 GeV. This 
requirement suppresses photoproduction background, in which the electron escapes detection 
and a hadron fakes the electron signature. Events are accepted if the z coordinate of the event 
vertex, reconstructed using the central tracker, lies within ±35 cm of the mean position for ep 
interactions. 

2.3 Selection of the Baryon Candidates 

The hypothetical doubly charged X and the neutral X° baryon states are identified by com- 
plete reconstruction of their respective decay chains through and A baryons into pions and 
protons, according to 

^Herein, the term "electron" is used generically to refer to both electrons and positrons. 
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X ^ S TT [AtT ]7r [{pTl )'K ]'K (1) 

X° ^ S-7r+ ^ [A7r-]7r+ ^ [(p7r-)7r-]7r+ . (2) 

The decay daughter particles are fitted in three dimensions to their respective decay vertices 
[fT9l . which are referred to as the tertiary (A-decay) and the secondary (H-decay) vertex. The 
analysis is based on charged particles measured in the central region of the HI detector with a 
minimal transverse momentum of 0. 12 GeV. These (non vertex-fitted) tracks may in addition 
be kinematically constrained to the primary vertex, in which case they are specifically referred 
to as primary vertex-fitted tracks. The baryon selection is chosen to optimise the signal for the 
standard baryon H(1530)° . 

In the first step, the A baryons are identified by their charged decay mode, A ^ pn^ , using 
pairs of oppositely charged tracks. The track with the higher momentum is assigned the proton 
mass. The particles are fitted to the tertiary vertex and the A candidates are retained if the 
vertex fit probability is above 1%. To reduce background, the two-dimensional separation (in 
the plane transverse to the beam line) of the tertiary and the primary vertex is required to be 
larger than 7.5 mm. The transverse momentum of the A is required to be larger than 0.3 GeV. 
The contamination from K'g vr+vr" decays is suppressed by excluding candidates with a tttt 
invariant mass in a ±10 MeV (~ 2a) window around the nominal Kg mass. A candidates are 
retained for the further steps if their invariant pn mass is within ±8 MeV of the nominal A mass, 
which is subsequently assigned to the selected A candidates. 

In the second step, candidates are formed by combining each of the A candidates with 
a negatively charged track assumed to be a pion. It is required that: the A candidates and the 
negatively charged tracks be fitted to a secondary vertex with a probability above 0.1%; the two- 
dimensional distance of closest approach in the transverse plane (DC A) of the candidates 
with respect to the primary vertex be smaller than 2.5 mm; the angle between the secondary 
and the tertiary vertex vector^ be less than 0.6 rad. The invariant mass spectra Af(A7r^) and 
M(A7r+) of all candidates passing these criteria are shown in figure [IJ Clear signal peaks are 
observed around the nominal S" mass. S candidates are retained for the third step if their 
invariant mass is within ±15 MeV of the nominal mass, which is subsequently assigned to 
the selected candidates. 

In the third step, X baryon candidates are formed by combining each S~ candidate 
with one additional primary vertex-fitted track assumed to be a pion. For these combinations it 
is required that: the additional n"^ track have a significance (DCA/aDCA) to the primary vertex 
of less than 4.0; the transverse momentum of the X /° candidates be larger than 1.0 GeV and 
their pseudorapitity be within |?7| < 1.5. 

2.4 Simulation of Baryonic States 

To estimate the acceptance, the efficiency and the resolution for the detection of a hypothetical 
baryon state, a Monte Carlo simulation based on the PYTHIA [20| event generator is used, 

vertex vector is defined as the vector pointing from the primary vertex to the given vertex. 
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incorporating the Lund string model of fragmentation [llTTl . The kinematic distributions of 
strange S=—l hadrons in DIS data have been found to be reasonably well described 112211 by the 
PYTHIA simulation. A generic X baryon state is introduced in the simulation by changing 
the mass of the known A antibaryon to values in the required range from 1600 — 2300 MeV 
and having it decay into S^vr^ final states. Because most theoretical predictions and the mea- 
surements [|4l suggest a very small intrinsic width cr < 1 MeV for pentaquark states, the width 
for the X is assumed to be zero in the simulation. Due to lack of knowledge of the pentaquark 
production mechanism, it is assumed that the production kinematics of the X are similar to 
those of the standard H( 1530)° baryon. The generic X° state is introduced by changing the mass 
of the known S(1530)° baryon to values in the required range from 1600 — 2300 MeV, and let- 
ting it decay into S^7r+. All the generated events are passed through a detailed simulation of 
the HI detector response based on the GEANT program [|23ll and the same reconstruction and 
analysis algorithms as used for the data. 



3 Mass Spectra and Limit Calculation 

The resulting invariant mass spectra for each of the two neutral (S^7r+, S+tt^) and the two 
doubly charged combinations (S^tt^, S+7r+) are shown separately in figure [2l In both neutral 
spectra the signal of the well known S(1530)° state is observed. The sum of the two neutral 
and of the two doubly charged mass spectra are shown in the upper part of figures 3 and 4, 
respectively. 

The ratio of the invariant mass spectra of the neutral to the charged states is consistent with 
being flat in the mass region above the S(1530)° mass. This justifies the use of the same back- 
ground shape function for both the doubly charged and the neutral combinations. A simultane- 
ous fit of the neutral and the doubly charged mass spectra, shown in figures |3]and|4l is performed 
using a function F, that contains a Gaussian function G for the signal S( 1530)° baryon and a 
function B for the background shape, according to 



F = G + {l + Po)B; B{M) = Pi{M -m^-m^f' X (1 + PsM + PiM^). (3) 

Here, M denotes the Htt invariant mass and iris and m^r the masses of the H and the tt, respec- 
tively. The normalisation, the central value and the width of the Gaussian function G and the 
parameters P^ are left free in the fit. Pq represents the relative normalisation of the neutral to 
the doubly charged combinations and is set to zero for the neutral combinations. The fit yields 
a total of 163 ± 24 (stat.) 2(1530)° baryons. The reconstructed mass of 1532 ± 2 (stat.) MeV 
is consistent with the nominal value lfl?ll . The measured width of 9.4 ± 1.5 (stat.) MeV is in 
agreement with the detector simulation. 

No signal of a new baryonic state is observed above the H(1530)° mass in either the neutral 
or the doubly charged mass spectra (figures |3]and|4l). The resonance search can also be per- 
formed relative to the observed signal of the known S( 1530)° baryon, using the ratio R, which 
is defined as 
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R(M) = ^XP^ X (4) 

^ ^ A^(1530,0) e(M,g) 

iV(1530, 0) represents the number of observed S(1530)° H^7r+ and their antiparticle decays. 
N'^^^{M,q) describes the estimated number of resonance decays depending on the mass M 
and the charge q of the final state, which is derived from the difference between the observed 
spectra and the expected background contribution. The background distribution is taken to be 
the fitted function given by equation [3] For the calculation of N"^^^, the mass distribution of 
the signal is assumed to be a Gaussian function with a mean M and a mass -dependent width 
cr(M) corresponding to the experimental mass resolution. This width cr(M) varies from 6.8 
to 22.8 MeV in the mass range considered here, as obtained from the Monte Carlo simulation. 
The term e(M, q) describes the reconstruction efficiency of the S tt final state as determined by 
Monte Carlo simulation, and it depends on the Stt invariant mass M and on the charge q of 
the final state. Accordingly, e(1530, 0) represents the reconstruction efficiency for the neutral 
S(1530)° baryon. The ratio of efficiencies in equation |4] compensates for the small difference 
in the reconstruction efficiencies of the S(1530)° baryon and of a hypothetical baryon state. 
This factor varies from 0.96 (0.92) at the smallest masses to 0.87 (0.86) at the highest masses 
considered for the neutral (doubly charged) states. 

The ratio R{M) has the advantage that the systematic effects of the acceptances and the 
reconstruction efficiencies mostly cancel, making it insensitive to detector effects. In the ab- 
sence of an observed signal, the ratio R{M) therefore provides a robust quantity to set upper 
limits on the production of new narrow baryonic resonances decaying to S^vr''^ in the mass 
range 1600 — 2300 MeV. A mass-dependent upper limit at the 95% confidence level (C.L.) on 
the ratio R{M) is obtained from the observed invariant mass spectra using a modified frequen- 
tist approach based on likelihood ratios ll24ll analogous to the one applied in [51 . The method 
takes the statistical and systematic uncertainties in the number of signal events and background 
combinations into account. 

The systematic uncertainties on the ratio R{M) include the following contributions: 



1. The uncertainty on the overall number of S(1530)° candidates is 15%, as determined 
from the fit to the mass spectra. 

2. The statistical uncertainty on the determination of the reconstruction efficiency e(M, q) by 
the Monte Carlo simulation is taken as the value of the systematic uncertainty on the ratio 
e(1530, 0)/e{M, q). It is found to be below 8 % for all masses and charge combinations 
considered. Systematic uncertainties due to the lack of knowledge of the production 
mechanism of the hypothetical new baryon state are not taken into account. 

3. An uncertainty on the experimental width of the potential pentaquark signal is estimated 
to be 5 %, as determined by the difference of the observed width of the S(1530)° signal 
in data and simulation. 

4. A systematic uncertainty on the background distribution is assessed by performing the fit 
under different assumptions: using the background function ^ in the full mass range, 
excluding a mass window of ±2a around the considered mass, and also using the sum 
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of the background function and a Gaussian with fixed mass M and width a{M) to ac- 
count for a possible signal. The uncertainty of the number of background combinations is 
estimated from the difference between the different fitting methods and amounts to 2 %. 

In the limit calculation, the contributions from items 1-3 are quadratically added and included 
in the uncertainty on the signal events, while the contribution from item 4 enters the background 
uncertainty. 

The final results of the limit calculations are quoted for the kinematic region 2 < < 
100 GeV^ and 0.05 < y < 0.7, forp^(S7r) > 1 GeV and |r/(H7r)| < 1.5. As stated before, 
it is assumed that new resonances are produced by a similar mechanism to that of the known 
J = 3/2 baryons, that they decay into Svr with a 100% branching ratio, and that their natural 
widths are below the experimental resolution. 

In figure[3]the S tt^ invariant mass spectrum summed over the two doubly charged combi- 
nations is shown. In the same figure the 95% C.L. upper limit on the ratio R{M) is presented 
as a function of the mass M. The non-observation of a resonance state in the mass range 
1600 — 2300 MeV limits the production rate of a hypothetical "E^^ pentaquark to 12 — 45% of 
the S(1530)° production rate at the 95% C.L., depending on the (S7r)-mass. 

Furthermore, no signal is observed in the neutral invariant mass spectrum (figure ^ in the 
mass range 1600 — 2300 MeV, above the S(1530)°baryon. Figure |4] also shows the resulting 
upper limit on the ratio R{M), which is obtained in a procedure analogous to the doubly charged 
case. This restricts the production rate of a hypothetical pentaquark state to less than 10 — 
50% of that of the S(1530)°baryon, depending on the (S7r)-mass. 

This result is similar to the 10 — 50% upper limits measured by the ZEUS Collaboration ifTll 
in the region > 1 GeV^. 



4 Conclusion 

A search for new narrow baryonic resonances decaying into S^vr^ and S^7r+ and their charge 
conjugate states is performed with the HI detector using a DIS data sample in the kinematic 
region 2 < < 100 GeV^ and 0.05 < y < 0.7, corresponding to a total integrated luminosity 
of 101 pb^^. The established 2(1530)° baryon state is observed. No signal of a new baryonic 
state is found in the mass range 1600 — 2300 MeV. Therefore mass dependent upper limits at 
the 95% C.L. are set on the production ratio of hypothetical states, such as E'^^ and 2°^ , to the 
total number of observed S(1530)° baryons. 

The results reported here are similar to the limits measured by the ZEUS Collaboration. 
The overall HI statistics in the samples is comparable with the NA49 collaboration's data. 
The HI limits obtained at HERA do not confirm the NA49 observation of potential pentaquark 
states. 
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M(A 71) [GeV] M(A n*) [GeV] 



Figure 1: The invariant mass spectra for a) Att^ and b) Att^ particle combinations. The solid 
lines show the result of a fit to the data using a Gaussian function for the S~ signal and a 
background function as defined in equation [3] (with replaced accordingly by ttia), while the 
dashed lines indicate the background function only. 
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Figure 2: The invariant mass spectra of the En particle combinations for the charge combina- 
tions a) S^7r+, b) S+7r^, c) S^tt^ and d) S+7r+. For the purpose of illustration, the solid lines in 
a) and b) show the result of a fit to the data using a Gaussian function for the H( 1530)° signal and 
a background function (dashed lines) as defined in equation[3l The same background parametri- 
sation is used in c) and d). 
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\ H1 Data 



2<Q==<100 GeV^ 




M(5"7i",5V) [GeV] 



Figure 3: The invariant mass spectrum for the doubly charged combinations S^vr^ and S+7r+ 
(upper part). The solid line shows the result of a fit to the data using a background function as 
defined in equation [3l The lower part shows the 95% C.L. upper limit on the ratio R{M) as a 
function of the mass Af , as defined in equation |4l 
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Figure 4: The invariant mass spectrum for the neutral combinations S~7r+ and S+vr" (upper 
part). The solid line shows the result of a fit to the data using a Gaussian function for the 
H(1530)° signal and a background function (dashed line) as defined in equation [3l The lower 
part shows the 95% C.L. upper limit on the ratio R{M) as a function of the mass M, as defined 
in equation m 
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